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Abstract. 

Employing a relativistic version of a hypervirial result, recurrence relations for 
arbitrary non-diagonal radial hydrogenic matrix elements have recently been ob- 
tained in Dirac relativistic quantum mechanics. In this contribution honoring Pro- 
fessor Lowdin, we report on a new relation we have recently discovered between the 
matrix elements (2|r^|l) and (2|/3r^|l) — where /? is a Dirac matrix and the numbers 
distiguish between different radial eigenstates — that allow for a simplification and 
hence for a more convenient way of expressing the recurrence relations. We addition- 
ally derive another relation that can be employed for simplifying two center matrix 
element calculations in relativistic atomic or molecular calculations. 
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Introduction 

The evaluation of expectation values is always required for relating quantum 
calculations to experimental results in atomic or molecular physics. In most cases 
such expectation values can be expressed as matrix elements of powers of a radial 
coordinate r; this comes about since these powers can be regarded as either exact 
terms in a certain potential (as in the Lennard-Jones, the Casimir or the London 
potentials) or as terms in a multipolar expansion of the interaction with the electro- 
magnetic field [1-5]. Matrix elements of that sort can be also regarded as starting 
points of certain useful approximation schemes, like variational or Hartree-Fock or 
configuration interaction methods, to which the late Professor Lowdin made im- 
portant contributions [6-14] and which admit relativistic extensions [15-17]. The 
treatment of electromagnetic interactions in the realm of multiphoton transitions 
in very intense laser fields usually needs a full quantum electrodynamics treatment 
[1], but this is time consuming. As a good approximation we can use instead the 
relativistic Dirac quantum mechanical formalism [2] . One has only to remember the 
successes this theory has achieved in dealing with the hydrogen atom [18,19]. One 
can start with the known states of the relativistic hydrogen atom [18,20-23] and then 
proceed, as in non-relativistic quantum mechanics, to expand the states of interest 
in terms of the former [24] . The problem is thus reduced to the evaluation of matrix 
elements of powers of r between relativistic eigenstates of the hydrogen atom — a 
much simpler problem. The bad news are that these evaluations also become cum- 
bersome. This calls for techniques adroit for evaluating them. In nonrelativistic 
quantum mechanics these techniques in the form of algebraic methods, recurrence 
relations or clever uses of hypervirial theorems abound (see, for example, [25-31], 
but in the relativistic domain they are rather scarse (but see [4,32-34]). Such lack 
of techniques is also manifest in the non existence of relations valid for two center 
matrix elements. 

The matrix elements of different powers of r between Dirac eigenstates of the 
hydrogen atom referred to above require the evaluation the following type of integrals 



where /(r) is any function of r, the kets \nje) stand for a bispinorial Dirac radial 
eigenstate of the hydrogen atom and the subscripts label different states. Any of 
such states, when projected on the |r) basis, become [2] 




(1) 
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1 / Fnjeir) \ 

{r\nje) = -\ . (2) 

\iGnje{r) J 

In equation (1) we used the shorthand Fk = Fn^jf.eki.f) and Gk — Gn,,j,,€k{'>^)^ that 
are called, respectively, the big and small components of the bispinor (2), and that 
are the solutions of the radial Dirac equation. This radial equation can be written 
[35] as (iJfc — Ek)il)k(r) = 0, where the are the energy eigenvalues, 

Hk = car[pr - ipekUk + l/2)/r] + ^c^ + Vk{r), (3) 

is the radial Dirac Hamiltonian, k is just a label — useful in what follows — , and 
f3 = diag(l, —1) is a Dirac matrix (for a pedagogical discussion of the radial equation 
and of the hydrogen atom in Dirac quantum mechanics see [36]), and [37-38] 

Writing the radial equation in a more explicit form, we have 

c2 + Vk (r) - Ek ctk Uk + 1/2) /r - cd/dr 

_ cek {jk + 1/2) /r + cd/dr -c^ + Vk (r) - Ek 

where n = 0, 1, 2, . . . is the principal quantum number, j = 1/2, 3/2, 5/2, . . . the 
total (orbital plus spin) angular momentum quantum number, e = (—1)-^+'^^/^, 
I [= j it 1/2, according to whether / refers to the big or to the small component 
of the hydrogenic spinor) is the orbital angular momentum quantum number, and 
Vfc(r) is any radial potential (of scalar type). Please note that here as in all of the 
paper we are using atomic units: nie = e = h = 1. The quantum number e equals 
+1 when I = j + 1/2 and equals —1 when I = j — 1/2 and it is related to the 
often used eigenvalue, k, of the operator /5(1 + S • L) by «; = —e{j + 1/2), where 
E = diag(o", <r) and a = {cTx,o'y,az) is the usual 3-vector spin operator. Notice also 
that, as in equations (2), (4), and (5) above, we can work in the somewhat easier 
to handle 2 dimensional subspace of the 4-dimensional Dirac operators. This choice 
has no relevance whatsoever for the final results. Closed forms for the integrals (1) 
or certain approximations thereof have been given in the Appendix of [38] and in 
[34,39,40] but even these become very cumbersome for many uses. 

On trying to overcome the complications mentioned above, recurrence rela- 
tions which can be used to compute general — not necessarily diagonal — matrix 
elements between relativistic hydrogenic states of have recently been obtained 



0, 



(5) 
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[37,38]. Given such recursions, it is only needed to evaluate at most 6 matrix ele- 
ments, 3 for and 3 more for (3r^, for obtaining every other matrix element of 
or of (3r^ between hydrogenic states. 

In this contribution honoring the memory of Professor Lowdin and his work, we 
want to review these recurrence relations and to discuss one more relation, previously 
unnoticed to us, which can greatly simplify the use of the recurrence relations already 
reported. Other purpose of this contribution is to derive an hypervirial-like formula 
that can be useful for evaluating relativistic two center matrix elements. 

The previously known recurrence relations 

The recurrence formulae have been obtained rederiving in the relativistic realm 
a non-relativistic hypervirial result previously used to derive the Blanchard recur- 
sion relations between hydrogenic matrix elements of r" terms [41,42]. The non- 
relativistic Blanchard relation has been found so interesting that it has been gener- 
alized [39]. We have to point out that our relativistic approach is totally different 
from Blanchard's; we start constructing an hypervirial and then proceed with several 
relativistic identities that have led us to the recurrence relations we are after [37,38]. 
This approach is inspired in a technique previously employed for similar purposes in 
non-relativistic quantum mechanics [41]. 

To obtain a relativistically valid recurrence relation, let us first compute matrix 
elements of the radial function ^(r) = H2f{r) — /(r)iJi, where Hi is one of the 
Hamiltonians — these are really the same Hamiltonian just evaluated in any of the 
two states 1 or 2 — appearing in Eq. (3), to obtain. 



E (n2j2e2|/(r)|nijiei) = {n2j2e2\H2f{r) - f{r)Hi\ni jiCi) 

= -ic{n2 j2^2\(y.r (j'(r) + ^^f^fif)^ \nijiei), ^ 

where the primes are used to indicate r-derivatives and we have introduced the 
symbols A^^ = £2(2^2 + 1) ± ei(2ji + 1) and = E2±Ei, —in Eq. (4) only A21 
and E~ are used, Aj^ and will be used later on. Then computing the matrix 
element of H2^{r) - ^{r)Hi, we get [37,38] 



{E r(n2i2e2|/(r)|ni jiei) = 
2r 



This formula is the relativistic extension of the non relativistic second hypervirial 
introduced in a previous work [41]. Notice that though the non relativistic hypervirial 
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sufficed in determining the non relativistic recurrence relation it is not so in the 
relativistic case, we will also need the following two results: 
First, 



E+E (n2j2e2|r^|ni ji ei) = (n2j2e2|c^ 



4 2 

2Z [icarr^-^{l - A) - E'r^-^] - E+ Xicarr^-'^lm ji ei). 



(1 - X)(3 



„A-2 



(8) 



were, for obtaining this relation we essentially repeat the steps leading to (7) except- 
ing that, at the end, we evaluate the result on H2^{r) + ^{r)Hi. 
Second, we will also need 



E- 



'21 



{n2 32 ^2\{-ioLrr^ ^)\nijiei) = 



{n2 32 €2 1 - c(A - l)r 



A-2 



^21 ^ 



(9) 



the detailed steps for obtaining (9) can be found in [38]. In Eqs. (8) and (9) and in 
most of what follows, we have explicitly used /(r) = and V{r) — —Z/r. 

The first of the recurrence relations we are after follows from eliminating the 
two terms, 2iE'^ Xcarr^~^ and 2icarr'^~'^, from (7) and (8). In this way we may get 



Co(n2i2e2|r'^|niii ei) = ^ Q (n2 ^2 e2|r'^ ^\ni ji ei) +'^di {n2 j2 e2\/3r^ Vu'iei), 

(10) 



i=l 



i=2 



where the numbers q, z = 0, . . . 3, are given by 

E+{E-)^A^, 



Co 

Cl 

C2 = C 
C3 = 



D 



2Z{E-yA 



-^2 A- 
21 



^ A^iA^i c^A2-iA(A-l)E+ 

4 D 
2Zc2(A- l)(A-2)A^i 



L> + 4c2 



and the numbers di, i = 2 and 3, by 



(1 - A) + 



D 



d3 = 



Zc^(A- 1 )A 

Q' 



21 



(11) 



(12) 



where, for the sake of conciseness, we introduced the symbols 

/^2iE^ -4c2(A- 1) 



D = A^^E- 



4c2A, 



Q = 



(13) 



The second recurrence relation can be obtained from (8) through the simple but 
somewhat contrived process explained in [38], that yields 



eo{n2 j2 e2\(3r^\ni ji ei) = bo{n2 j2 e2\r^\ni ji ei) + 62(^^2 ^2 £2 ^|?T-iiiei) 
+ ei(n2 j2 e2|/3r^~^|ni ji ei) + 62(712 ^2 e2\l3r^~^\ni ji ei), 

where the numbers br and e,., r = 0, 2, 3, are given by 

6o=4A [{E-f-4c'] , 

b2 =c2(l-A) [(A2l)2-4A2], 

Co =2E+D, 

ei = - 4ZL>, 

e2 =c2^[(A,l)2-4A2]. 



(14) 



(15) 



Of course, the energy eigenvalues are those of the relativistic hydrogen atom, namely 



-1/2 



E = E 



= C 



1 + 



\ (ua - ja - 1/2 + V(ia + 1/2)^ - Z^aj.) J 



(16) 



where ap = 1/c ~ 1/137 is the fine structure constant. Please note that contrary 
to what it was implied in [37,38] the recurrence relations (8) and (12) remain valid 
even in the limit when A~ vanishes. 

The new recurrence relations 

Relations (9) and (12) can be useful for different computations, but they could 
be even more so if we could disentangle the matrix elements of from the matrix 
elements of ^r^. This can be achieved with the help of the following relationship 



At 

{n2 32 (^2\-^f{r)\niii d) 



2r 
{n2 32 €2! 



/3/(r)+^/'(r)|nijiei); 



(17) 



as follows from equation (8) by using the above illustrated technique but with dif- 
ferent substitutions. Notice that we have reverted to an arbitrary radial function 
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/(r) and a generic radial potential V{r). Now using, as we did before, the specific 
function /(r) = and the Coulomb potential V{r) = —Z/r, we obtain the useful 
new relationship 



go{n2 j2 e2\r^\ni ji ei) = jo{n2 j2 e2\Pr^\ni ji ei) + ji{n2 j2 e2\Pr^ Viii^i); (18) 
where 

go =A+ /2, 

JO =4^(^^i + ^2~i) + (A - 1), (19) 
h =^(A+ +A2-1). 

Equation (18) is fundamental for the disentanglement of the recurrence relations 
(10) and (14). 

Let us explain now the process for disentangling the matrix elements of 
from those of /?r^. We first obtain (71,2^2 ^2\i^^\ni ji Ci) from equation (16) and from 
equation (14); by equating the resulting expressions, we get the recurrence relation 
for terms of the form (712^2 ^2\Pr^\ni ji ei), as 

3 

?7o(«'2i2e2|/3r^|ni jiei) = ^ r?i(n2 ^2 £2! /3r^"' |ni ji ei), (20) 

i=l 

where the constants r]a, a = 0, 1, 2, 3, appearing in (19) are 



^0 



E+D 



R 



2A 



"2A((£-)2 -4c4) 2c2 A+' 



.R 



D 



V3 



(A-l) 4A^-(A,-,)^ 
2A (E-f-Ac^) 



E+R c2(A-2) c2A2i 



Z{X-1)R 
' 4A 



4A^ - i^2i? 
{E-f - 4c4 



A2+1 



4(A-1) 



(21) 



Equation (20) is, together with the definitions (21), the new recursion for matrix 
elements involving just /3r^. 

To obtain the recursion for the matrix elements of r", we start with equations 
(10) and (18). We first obtain the three (^2 j2 e2|/3r'^~*|ni ji ei) {i = 0,1,2) terms 
from them; next we succesively substitute them into equation (14) and, by juggling 
with the resulting equations, we are able to obtain 
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i^o(?T'2j2e2|r-^|ni ji ei) = ^ i^i (n2 j2 e2|r^ InijiCi) 



1=1 



where the constants appearing in (22) are 



z/i = - 8£;+(£; 



Q 



Zc2(A- 1)' 

_,2(Q + 6c2) 



2Ar 



c2(A-l)' 



E+_ 
Z 



D 



J 



U3=-D 



8Z'^{E-fW+—E+{E-) 



2 4A2-(A2-i) 



2" ~ "r" 

21 



4£;+(A-2)gVF 



QW 

{E-f (4A2 - (A,,)2) 
(A-1) 



c\\-l) (4A2 - (A,-,)2) 



2Z 



8Z(A - 2)QW 



XE+ 



^ ^ 4c^(A-2) 
L(A-1)W^ ^ 



2Z 



(4A^ - (A,,)2) 



[/ _ c2 
4 R 



(22) 



(23) 



i.5=c^(A-2) (4A2 - (A^D^) . 
Again for the sake of conciseness, we have introduced the following definitions 



Q 



A^i^;- -4c2(A- 1) 



A2~l 



T = {{E-r-4c^), 



W = 



A^,E- 



4c2A 



U = 



Asi^;- -4c2(A- 1)' 
{A,-,^- - 4c2(A - 1)}[4A(A - 1)E+ - A+ (Aj^E- - 4c2A)] 
(A2-iE--4c2A)(A-l)A+ 



(24) 



Equation (22) is the recursion involving only matrix elements of r". 

The new recurrence relations are equations (20) and (22), we also think that 
equation (19) can be found useful in certain applications. We must also pinpoint that 
the relationships derived are valid in general for complex values of the exponent A, as 
long as cJi +0^2 + 1 A I > — 1 , were the numbers cui are defined as + 1/2)2 _ ^2^2^_ 
This property was established in [38]. 

A relation for two-center integrals. 

The recurrence relations established above for relativistic hydrogenic states is 
sufiiciently general for performing many calculations. It is not however useful for 
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calculations involving matrix elements of a radial function between states represented 
by radial eigenstates corresponding to different potential functions, that is, the so- 
called two center matrix elements. It is the purpose of this section to present one 
such relation for relativistic two-center matrix elements that can be regarded as a 
first step in the appropriate direction. This will be done using ideas developed in 
[38,41,43]. 

General hypervirial results, the virial theorem, and other algebraic techniques 
[25,26,44] have been always very useful for calculating matrix elements, nonetheless 
they have been little used in relativistic calculations (but see [37,44,45]). We want 
to show here that they can have comparable importance than in non-relativistic 
quantum calculations. So, let us consider two radial Dirac Hamiltonians [as in equa- 
tion (3)] with two possibly different radial scalar potentials Vi{r) and V2{r) — we are 
thinking on potentials like those describing vibrational states in a molecule — 

Hi=car[pr-ipei{ji + l/2)/r]+^c^ + Vi{r), 

H2 =car[Pr - iPe2{j2 + l/2)/r] + Pc" + V2(r), 

a further assumption we make is that the two potentials are displaced from their 
respective equilibrium points by a constant quantity, i. e. ri + r2 = a, where 
is the equilibrium position of potential Vi{r) and a the displacement. Notice that, 
at difference with equation (3), here the numerical label in the Hamiltonians is not 
just convenient but has a definite physical meaning, for H2 and Hi refer to two in 
principle different systems. 

Taking the difference between the Hamiltonians (25), Hi — H2, we get 

Hi=H2 + icar(3^ - {V2 (r) - Vi (r)) . (26) 
On employing (25), we can inmediately evaluate the commutator 

[Hij2{r)] = -icar^^ (27) 

where /2(?") is an arbitrary radial function and [H,f{r)] stands for the commuta- 
tor between H and /(r). We can calculate this commutator again, but now using 
expression (26) to get the alternative form 



[Hi, /2(r)] = H2f2ir) - f2{r)Hi - (F2 - Vi) /2(r) + ica,/3^/2(r). (28) 

If we now equal (28) with (27) and take matrix elements of the resulting ex- 
pression between the two states, (nijiCil and 177.2^2^2) — which correspond to the 
different Hamiltonians (25) — we directly obtain the hypervirial inspired relation 
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{E2 - El) {mil ei|/2(r)|n2 j2 £2) = (^1 Ji ei| (^2 - Vi) /2(r)| 712^2 £2) 
-zc(m jieila,. ^/2(r) +/3^/2(r)^ |n2i2e2>; ^ 

notice, however, that the energy eigenvalues in (29) refer to the different Hamiho- 
nians (26), and not to different states of the same Hamihonian — as was the case 
analysed in the first section of the paper [compare with equation (6)]. In fact, if we 
consider the same potentials in equation (29), i. e. we take V2{r) = Vi(r), we are 
also inmediately setting Ei = E2 and thus we recover the relationship (6), which is 
valid for just one center. 

The important point is that equation (29) is an exact relation valid for the 
calculation of two center matrix elements, given as a function of the eigenenergies 
for any two scalar radial potentials Vi{r) in the Dirac equation, and thus it can be 
useful for deriving recurrence relations between such matrix elements. But this is 
still a work in progress. 
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